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Abstract 

The work describes a possible physical characterization for the defi- 
nition of elementary fluid flow. As consequence, an analytical expansion 
based on Reynolds number powers for the velocity field is shown in the 
weakly turbulent case. 

1 Notion of elementary analytical flow 

Let us recall the notion of elementary analytical flow (see pQ). 
Let L a length unit and Cl a closed domain of M n , n > 2, formed by a small 
number (e.g. 2 or 3) of L 2 squares or L 3 cubes. In a 2D steady case, let F a 
fluid entering in £1 from a single edge of a square and flowing out from another 
single edge. Let also s, < s < 1, a parameter which describes the position of 
a single fluid particle along the inflow edge, so that sqL is the initial position of 
the generic particle associated to the particular value sq of the parameter. 
Let s i — > R(s) a function R : [0, 1] — > [0, 1] which maps a value of the pa- 
rameter s to the value identifying the position reached by the particle on the 
outflow edge, so that this position is identified by the value R(s)L. If x, y is 
a cartesian coordinates system in the plane of £1, assume that each streamline, 
or each particle path in lagrangian view, is described by a parametric curve 
t i ► ^(t) = (x(t), y(r)), with < r < 1. The parameter r could not be, in 
general, the time variable of the flow. Let this parametric representation be 
determined by a set A of analytical conditions regarding $(r) and $(r), that 
is the passage of the particle in some suitable points of 51 and the velocity field 
direction in some other (or the same) points. 

An elementary analytical flow is then a particular list F = {O, R, A}. 



In this work let f2 the domain formed by three L 2 squares, the first two along 
the x-axis of a cartesian coordinates system from to 2L, and the third above 
the second, from L to 2L y-coordinates (PJ)- 
Now the set A af analytical conditions is so defined: 

Pi. at the inflow edge, for r = and position parameter s, let $(0) = (0, sL); 

Pi- at the outflow edge, for r = 1 and position parameter R = R(s), let 
*(l) = (L(l+iZ),2L); 

P3. for r = i a particle path intersects the diagonal line of the second square, 
that is the line of cartesian equation y = —x + 2L; if a parameter p, < p < 1, 
describes the positions along this line, the condition is $(5) = (L(2 —p),pL); 

Di. at the inflow edge the velocity is parallel to x-axis, so that for every s 
y(0) = 0; 

£>2- at the outflow edge the velocity is parallel to y-axis, so that for every s 
x(l) = 0. 

Assuming that the components of $(t) are cubic polinomials on variable r, that 
is the more simple analytical expressions satisfying the set of conditions A, by 
usual computations follows (PP) 



x{t) = L(10 -6R- 8p)r 3 + L(-21 + 1LR + 16 P )t 2 + L(12 -4R- 8p)r (1) 



y(r) = L(4 + 6s 



)r 3 + L(-2 - 7s + 8p)r 2 + Ls (2) 



Also, let p — s, so that particles paths have their initial shape until the line 
y = —x + 2L. In this case the curve $ has the form 



x(t) = L(-6R -8s+ 10)r 3 + L(UR + 16s - 21)r 2 + L(-AR - 8s + 12) (3) 

y(r) = L(2 - s)(2t 3 - t 2 ) + Ls (4) 

Assume now the following definition for R(s), which gives a possible mathemat- 
ical modelization of the physical phenomenon of weak turbulence (see pQ and, 

e.g., m) 



1 + sin I 1 Me s 

cw 



(5) 



where Me is the Reynolds number and cp is a constant associated to the flow. 
The previous form of R(s) is a consequence of the following Taylor expansion 
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sm(l - 2s)] « 1-6 
which is the correct shape of R(s) in the laminar case ([I])- 



lim R(s) = - [1 

£e — ►cf & 



(6) 



2 A possible expression for cjf 

The relation is the viscous limit (the inverse of the inviscid limit, see (0)), 
which is valid for small values of Reynolds number. On the contrary, when Re 



2 



is large, the argument of the sine function in (J5J should assume great values 
for the frequency so that the flow becomes turbulent because, for a fixed 
particle, R(s) can be a value very different compared to laminar value 1 — s. 
If we would fixed the constant cf indipendently from the flow considered, that 
is from the analytical representation $(r), the most simple choice which satisfy 
the two previous asymptotic constraints is 

Substituting (JJJ into (JSJ, the viscous limit © is expressed by the relation 

lim Re 2 = 1 (8) 

Re-*c F 

that is Cf = Me = 1, which can be viewed as a physical definition of viscous 

limit for elementary flows. 

A more general choice could be 

cw = -=rw— r ( 9 ) 

but in this case we should estimate the coefficients Ai , and the value of Re in 
the case of viscous limit could be not unique. 



3 Choice of and time parameterization 

In this section we show a relation between cf and time parameterization for an 
elementary flow. 

Assume that the parameter r of the paths representation <5(t) could depend, 
at least as first approximation, linearly on time variable t by a multiplicative 
regular function / = /(s), where s is the parameter identifying a single particle 
path from inflow position: 



T = f( S )t (10) 

In this way, the flow velocity field (ui(s, t), U2(s, t)) can change from one path 
to another. Note that, for a fixed value of s, this formula has meaning only for 
< t < f(s) . For every fixed s follows 



dx dx dr , , . 

u 2 = — = — — = f(a)y 
dt dr dt 

Now, let t — 0. Then, being in the considered example y(Q) — Vs, the 
initial velocity field at inflow edge, usually a known boundary condition, is 
(ui(s,0),n 2 (s,0)) = («i(s,0),0). From © follows that i(0) = 4£(3 - R - 2s), 
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and at least for those inflow points for which R(s) ^ 3 — 2s the analytical 
expression of f(s) is 

The analytical relation between velocity field and particles paths is then 

( Ul ( S ,t),u 2 ( s ,t)) = (i(r),y(r)) ^ t (13) 



4 An expansion for velocity field 

Using (JJJ, suppose now that Reynolds number Re is small enough to have 
1 - 2Re 2 s w 0, so that at first order R{s) » 1 - Ke 2 s. Then x(0) = 4L[2 + 
(Re — 2)s]. Also, assuming regular properties for Ui(s, 0), from Taylor expansion 
centered e.g. on s = let Mi(s,0) «io + uiqs + ^ilios 2 , where Uxo = ui(0, 0), 
"io = d s ui(0,0) and u\o = <9 2 s wi(0, 0). Using a Taylor expansion for /(s), it 
can be shown that 



/(*) 
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(14) 
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This formula shows the role of Reynolds number and inflow boundary condi- 
tion about the time parameterization l|10(l of the particles paths and about the 
expression of velocity field (0U+G3)- For high values of Re, the expansion 
of -R(s) should consider higher order terms, so that the approximated expres- 
sion for the velocity field could present higher power of Reynolds number (for 
a method of expansion based on Re powers see (0). Usual computations give 
the following series expansion for /(s), in the case of a parabolic inflow velocity 
«i(s, 0) = v + V\s + v 2 s 2 - 
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(15) 
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Using the two relations 1)11(1 . previous formula gives an expansion of velocity 
field in terms of Reynolds number powers. In the attached figures we present 
some profiles for the scalar speed \/ u\ + u 2 in the case of an initial Poiseuille 
flow ui(s, 0) = 5(s — s 2 ). 
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T = 0. 
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